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1 Introduction and statement of result 

PT-symmetric operators in quantum mechanics (for general references on VT- 
symmetry and quantum mechanics, see e.g. [3]) are those operators on R" left 
invariant by a successive application of the parity operator P, acting on wave func- 
tions as 

{ViIj){x) = ip{{-iy^xi, . . . , {-ly^Xn), 3k = 0, 1, not aU equal to 0, 



and of the time-reversal symmetry, acting as {Ti/j){x) = ip{x). The PT-symmetry 
of an operator is called exact if its spectrum is purely real, see e.g. [1]. Generally 
speaking, the reality of the spectrum, and thus the exact "PT-symmetry, cannot 
follow by a unitary equivalence to a self-adjoint operator, since a "PT-symmetric 
operator is not necessarily self-adjoint, and actually, the standard examples are not 
self-adjoint and not even normal operators. However, it could follow by a similarity 
to a self-adjoint operator. Let us therefore proceed now to discuss the notion of 
similarity for two unbounded linear operators. For the limited purposes of this 
paper, here we shall only consider a rather particular abstract situation, introducing 
assumptions which will be satisfied in the specific instances below. 

Let T-Lj, j = 1,2, be complex separable Hilbert spaces, and let Aj : Tij — Tij be 
closed densely defined operators such that for j = 1,2, the spectrum Spec(^j) C C 
is discrete, consisting of eigenvalues of finite algebraic multiplicity. When A G 
Spec(^j), we let e{^^ C ^{Aj) be the finite-dimensional space of generalized eigen- 
vectors, corresponding to Aj, A, so that 

E?) = Ker (A, - Xf , 
if > N{X,j) is large enough. 
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Let Sj C V{Aj), j = 1, 2, be linear subspaces such that Aj{Sj) C Sj, and E^^' c Sj, 
for each A G Spec(^j). Assume that S : Si ^ S2 is a hnear bijection such that 

SAi = A2S on Su (1.1) 

or equivalently, such that 

AiS-^ = S-^A2 on ^2. (1.2) 

We shall then say that the operators Ai and A2 are related by the similarity trans- 
formation S, or are similar. 

Proposition 1.1 The similar operators Ai and A2 are isospectral. 

Proof: Let u G E^^\ so that (^1 - X)^u = 0, for some A'^ G N. It follows that 
u, {Ai — X)u, ... E Si, and we can apply fll.ip to see that {A2 — X)'^Su = 0, so 
Su G Thus we have shown that S{e{^^) C E^^\ The same argument applied 

to S^^ gives E^^^ D S~\Ef^), and we get Spec(^i) = Spec(^2), and S{e[^^) = Ef\ 
for all A G Spec(^i) = Spec(A). □ 

Notice that in the discussion above, the operator A2 may be self-adjoint on 7^2, 
even though the operator Ai need not be self-adjoint on "Hi. It is therefore natu- 
ral to ask whether a PT-symmetric operator with real spectrum is always similar 
to a self-adjoint one. This problem has attracted considerable attention, notably 
by [22], [23], [23], in an abstract framework. Here instead, we shall consider the 
class of PT-symmetric elliptic quadratic differential operators acting on L^(R"), 
with real spectrum, and establish necessary and sufficient conditions for their simi- 
larity to a self-adjoint operator, as above. A noticeable peculiarity of these conditions 
is their classical nature, i.e. their dependence only on the classical flow generated 
by the classical quadratic hamiltonian, whose quantization yields the given VT- 
symmetric quadratic differential operator. Let us also remark here that in general, 
the operator S in (11.11) realizing the similarity between the non-self-adjoint opera- 
tor Ai and the self-adjoint operator A21 cannot be expected to be bounded with a 
bounded inverse, due to the issue of pseudospectra for non-self-adjoint operators, 
see [TO], [9], [2]. 

Let q{x,^) be a complex- valued quadratic form on R^" = R" x R". Here (x, ^) G R^"' 
are the canonical coordinates, so that {^i,Xj} = 6ij, I < i, j < n. Throughout this 
work, we shall assume that the quadratic form q is elliptic on R^", in the sense that 
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q{X) = 0, X G R^", if and only if X = 0. An application of Lemma 3.1 of [25] then 
shows, if n > 1, that there exists z G C\{0} such that Re (zq) is positive definite. In 
the case when n = 1, the same conclusion holds, provided that the range of q on 
is not all of C, which will be assumed in the sequel. After a multiplication by z, we 
shall assume in what follows, as we may, that 2; = 1, so that Reg is positive definite. 
It follows that the range E(g) := g(R^") of q on R^" is a closed angular sector with 
a vertex at zero, contained in the union of {0} and the open right half-plane. 

We shall now introduce the assumption of the "PT-symmetry of the quadratic symbol 
q. To that end, let n : R" — )■ R" be linear and such that 

= 1. (1.3) 
Associated to the involution k, we have the parity operator V, given by 

{Vu){x) =u{k{x)), (1.4) 
and the lift of k, to the cotangent space, given by the linear involution, 

;C:R^"->R2", IC{x,0 = {<x),-^\0)- (1-5) 
The PT-symmetry assumption on the quadratic symbol q is of the following form. 



qolC = q^ q{K{x), -K^{0) = q{x, 0, i^, ^ R'"- (1-6) 

It follows, in particular, that the sector S(g) is symmetric with respect to reflection 
in the real axis. 

Associated to the symbol q is the corresponding operator Op"'(g) acting on L^(R"'), 
formally defined as the Weyl quantization of q{x,^), 

Op^{q)u{x) = (27r)-" / / e*«-(--^)g (^,^) ^(y) dyd^ (1.7) 

Writing 

q{x,0= Yl (1-8) 

\a+f)\=2 

we see that the operator Op'"(g) is a differential operator of the form 

Op-(g)= Yl 

\a+l3\=2 
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Let us recall from that the maximal closed realization of the operator Op"'(g), 
i.e. the operator on L^(R") equipped with the domain 

{u e L\W)- Op'"{q)u e L2(R")} = {u e ^^(R"); x^Dfw G L2(R"), |a + /3| < 2}, 

agrees with the graph closure of its restriction to iS(R'*), 

Op"(g) : 5(R") ^ 5(R"). 

Furthermore, the spectrum of Op"'(g) is discrete, and its precise description will be 
recalled below. 

The assumption (11.61) implies the PT-symmetry of the operator Op"'(g), i.e. the 
commutation property [Op"'(g), P7~] = 0. Here the parity operator V has been 
introduced in (II. 4p . Indeed, we have 

V o Op'"(g) o = Op'"(g(«:(x), «:*(0)); T o Op'"(g) o = Op'"(g(x, -0), 
Vr o Op'"(g) o (VT)-' = Op'"(g(«:(x), -«:*(0))- 

It follows that the spectrum of Op"'(g) is symmetric with respect to the real axis. 

The essential role in what follows will be played by the fundamental matrix F of the 
quadratic form q. When recalling the definition of F following Section 21.5 of [18], 
we let 

be the canonical symplectic form on R^", which extends to the complex symplectic 
form on C^". Letting also q{X, Y) stand for the polarization of q, viewed as a 
symmetric bilinear form on C^", we define the 2n x 2n fundamental matrix F by 
the identity 

q{X, Y) = a{X, FY), X,Y E C^". (1.9) 

We notice that the fundamental matrix F is skew-symmetric with respect to a, and 
furthermore, following [18], we see that in the canonical coordinates {x, it is given 
by 

F=l('b % )• (1-10) 

We can now state the main result of this work. 
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Theorem 1.2 Let q : x C be a quadratic form, such that Reg is positive 

definite. Let n : R" — > R" be linear and such that = 1. Assume the property of 
the VT -symmetry , 

q{K{x),-K^{0) = q{x, 0, {x, e R'", 

and that the spectrum of the quadratic operator Op^{q) is purely real. Then the 
operator Op"'(g) is similar to a self-adjoint operator, in the sense of the discussion 
in the beginning of this section, precisely when the fundamental matrix F of q has 
no Jordan blocks. 

Remarks. 

1. Let us consider the Hamihon vector field of q, 

H, = q[-d,~q',-d^. (1.11) 
According to fll.lOp . we have 

FY = \h,{Y). 

Therefore the possibihty of estabhshing a similarity between a PT-symmetric 
elliptic quadratic operator with purely real spectrum and a self-adjoint oper- 
ator involves only the underlying classical fiow, i.e. is determined by a purely 
classical condition. 

2. Examples of PT-symmetric non-diagonalizable hamiltonians with real spec- 
trum can be found in [7], as well as in [1], [30], [8]. See also the discussion in 
Section 4. 

3. In Theorem 11.21 we have only considered the case when the quadratic form q 
is elliptic on R^". It seems quite likely, however, that the result of Theorem 
ll.2l can be extended to a suitable class of partially elliptic quadratic operators, 
namely the one studied in [H] , [31] . We leave this extension open until needed. 

The plan of this note is as follows. In Section 2, we establish a symmetry pro- 
perty of the fundamental matrix F, as a consequence of the assumption of the 
T^T-symmetry of q. The proof of Theorem 11.21 is then carried out in Section 3, 
using the techniques of FBI-Bargmann transformations with quadratic phases [26] , 
which, in the quadratic case, allow us to construct the similarity transformation 5* 
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explicitly. Section 4 is devoted to the discussion of an example, due to [8], of an 
elliptic quadratic "PT-symmetric operator with real spectrum, for which the funda- 
mental matrix possesses Jordan blocks. The appendix A gives a concise exposition 
of aspects of the theory of positive Lagrangian planes, required in the proof of the 
main result, while in the appendix B, the general theory of quadratic Fourier integral 
operators in the complex domain, arising when quantizing complex linear canonical 
transformations, is developed. 

Acknowledgements. Johannes Sjostrand acknowledges support from the French 
ANR grant NOSEVOL: ANR 2011 BS01019 01. We would also like to thank the 
referees for very helpful remarks and suggestions. This led to a generalization of our 
main result. 

2 P7"— invariance condition for the fundamental 
matrix 

In this section, we let q : R^" — C be a quadratic form such that Reg is positive 
definite. It follows therefore from (11. 9p that the eigenvalues of the corresponding 
fundamental matrix F avoid the real axis, and in general we know from Section 21.5 
of [IB] that if A is an eigenvalue of F, then so is —A, and the algebraic multiplicities 
agree. Let us also recall from [IB] that the eigenvalues of F belong to the set 
iS(g) U —iT,{q), and let Ai, . . . , A„ be the eigenvalues of F, counted according to 
their multiphcity, such that Xj/i G T,{q), j = 1, . . . ,n. From [6], [25], we know that 
the spectrum of the operator Op'"(g) is given by the eigenvalues of the form 



Notice that Spec(Op"'(g)) C S(g). 

Let us recall the real linear involution )C, introduced in (11.51) . Writing /C = J o /C^, 
where J{x,^) = [x, — ^) and Kiiy^rf) = {n{y), {n*)~^{r])) is symplectic, we see, using 
also that the involution J is skew-symmetric with respect to the symplectic form a, 
that the map K, is antisymplectic in the sense that 




(2.1) 



a{ICX,)CY) 



a{X,Y), X,Ye'R- 



In particular. 



a{}CX,Y) 



a{lCX,KY) 



a{X,}CY). 
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We shall also let C : C^" C^" stand for the involution defined by the operation 
of complex conjugation. We have the following result. 

Proposition 2.1 Assume that the quadratic form q with Keq > is such that 
q{]C{x,$,)) = q{x,^), (x,^) G R^". Then the fundamental matrix F of q satisfies 

F = -ICCo F olCC. (2.2) 

Furthermore, the eigenvalues of F are symmetric with respect to the imaginary axis. 

Proof: Let us write, 

q{x, ^) = Ax-x + 2Bx ■ ^ + C{ ■ ^. 

Here A, 5, and C are complex n x n matrices, with A and C symmetric. The 
condition of the PT-symmetry, 

g(/€(x),-/€*(0) = g(x,0, (x,0 eR'", 

implies then that the matrices A, B, and C satisfy 

A = k^Ak, B = -kBk, C = kCkK (2.3) 

When establishing the property (12.21) . we introduce the symmetric bilinear form 
associated to g, given by 

g((x, i),{y,r^)) = Ax-y + Bx-ri + Byi + a-r,, (x, G C^^, (y, r/) G C^". 

Using (12.31) . we see that 



g(/C(x, 0, /C(i/, r/)) = g(C(x, 0, C(y, r/)), 
and therefore, from (II. 9p we get 



a{CX,FCY) = a{ICX,F]CY) = -a{X,ICFICY) = -a{CX,CICFICY). 

Here we have also used that the map K, is skew-symmetric with respect to a. The 
identity ( 12. 2p follows, and it only remains to check the symmetry of the eigenvalues 
of F with respect to the imaginary axis. To that end, assume that A G C is such 
that FX = XX, X e C^"". li Y = JCCX, then (F o /CC) F = XJCCY. Applying 
the antilinear map ICC and using (12. 2p . we conclude that FY = —XY. The proof is 
complete. □ 

Remark 2.2. Combining Proposition 2.1 with the explicit description of the spec- 
trum of the PT-symmetric operator Op"'(g), given by ( 12. ip . we conclude that if the 
spectrum is purely real, then the eigenvalues A-,-, 1 < j < n, are purely imaginary. 
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3 Similarity transformation and proof of Theo- 
rem 11.21 

Assume that q : R^" — )■ C is a quadratic form such that Reg > 0, and let F be the 
fundamental matrix of g, introduced in fll.91) . When A G Spec(F), we let 

= Ker((F - A)2") c C (3.1) 

be the spectral subspace corresponding to the eigenvalue A. According to Lemma 
21.5.2 of [IB], the complex symplectic form a is nondegenerate viewed as a bilinear 
form on Ex x E_\. 

Let us introduce the unstable linear manifold for the Hamilton flow of the quadratic 
form i~^q, given by 

A+ := C C". (3.2) 

Im A>0 

According to Proposition 3.3 of the complex Lagrangian plane A"*" is strictly 
positive in the sense that 

-ct(X,X)>0, 7^XgA+. (3.3) 
i 

We refer to Appendix A for a discussion of positivity conditions in C^". Defining 
also 

A- = i^A C C^", (3.4) 

Im A<0 

we know from the arguments of [25] that the complex Lagrangian plane A~ is strictly 
negative in the sense that 

-ct(X,X)<0, Ot^XgA-. (3.5) 

It also follows from (11. 9p that the quadratic form q vanishes when restricted to A''^. 

Our proof of Theorem 11.21 will proceed by exhibiting a complex linear canonical 
transformation which will reduce A"*" to {(x, ^) G C^"; ^ = 0} and A~ to {(x, ^) G 
C^"; X = 0}. We shall then be able to quantize the canonical transformation by 
means of an FBI-Bargmann transform [26], which will essentially provide us with 
the sought similarity operator. Let us notice here that the use of such canonical 
transformations has a long and rich tradition, see for instance [21], [12], [28]. The 
discussion below will follow Section 2 of [15] closely. 



9 



Although not necessary, it will be convenient to simplify q first by means of a suitable 
real linear canonical transformation. When doing so, we observe that an application 
of Example A. 6 in Appendix A shows that the negative Lagrangian A~ is of the 
form 

7] = A.y, y G C", 

where the complex symmetric n x n matrix is such that lmy4_. < 0. Here 
{y,rj) are the standard canonical coordinates on T*Iiy = x R^, extended to 
the complexification T*Cy. Using the real linear canonical transformation {y,r]) 
[y, rj — [Ke A^)y) , we reduce A~ to the form rj = ilmA^y, and by a diagonalization 
of lm^_, we obtain the standard form rj = —iy. After this real linear canonical 
transformation and the conjugation of the quadratic operator Op^{q) by means of 
the corresponding metaplectic operator, unitary on L^(R"), we may assume that 
A~ is of the form 

r] = -iy, yeC\ (3.6) 

We refer to the Appendix to Chapter 7 of [11] for the introduction to and basic 
properties of metaplectic operators. Notice also that according to Example A. 6, in 
the new real symplectic coordinates, extended to the complexification, the positive 
complex Lagrangian A"*" is of the form 

ri = A+y, ImA+>0. (3.7) 

Let 

B = B+ = {l-iA+)~^A+, (3.8) 
and notice that the matrix B is symmetric. The holomorphic quadratic form 

ip{x, y)='-{x- yf - ^Bx ■ x, (x, y) G C^", (3.9) 

satisfies Imipyy > and detip'^y ^ 0. It gives rise therefore to the FBI-Bargmann 
transformation, 

Tu{x) = C I e^^(^'J')n(y) dy, x G C", C > 0, (3.10) 

and following [25], [21], the operator T is to be viewed as a Fourier integral operator, 
with the associated complex linear canonical transformation of the form 

: C^" 3 {y, -^'{x, y)) ^ (x, ^'Ax, y)) G C^". (3.11) 
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Using (13. Sp . ( 13. 9p . we see that is given by 

Kt : (y, ri) t-^ (x, ^) = - z?7, 77 + iBr] - By), (3.12) 

and therefore, we have kt (A_) = {(x,^) G C^"; x = 0}, while kt(A"^) is given by 
the equation {(x,^ e C^";^ = 0}. 

We know from [29] that for a suitable choice of C > in ( I3.10p . the map T is unitary, 

T: L2(R'^) ->i/^,,(C"), (3.13) 

where 

i7$„(C") = Hol(C") nL2(C";e-2*oL(rfx)), 
and $0 is a strictly plurisubharmonic quadratic form on C", given by 

$o(x) = supyeRu {-\m.ip{x,y)) = ^ ((Imx)^ + Im {Bx ■ x)) . (3.14) 

From [29], we recall also that the canonical transformation kt in (13. lip maps R^'^ 
bijectively onto 

A*o = {(^,^^(^));^eC"}. (3.15) 

The positivity of {(x,^) G C^"; ^ = 0} with respect to A$q implies, in view of 
Proposition A. 9 in Appendix A, that the quadratic form $0 is strictly convex, so 
that 

^o{x)r^\x\\ xgC". (3.16) 



Example: The standard Bargmann transformation. Following [3], let us 
consider a complex integral transform of the form (I3.10p . where the phase function 
ip{x, y) is given by 



V9(x,?/) = i 



The corresponding canonical transformation, defined as in (I3.1ip . is then of the form 
while the associated quadratic weight function is given by $(x) = (1/2) |x|^. 
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Returning to the FBI-Bargmann transformation given by (13. 9p . (I3.10p . let us next 
recall the exact Egorov theorem, |29j . 

TOp"'(g)M = Op"'(g)TM, u G 5(R"), (3.17) 

where g is a quadratic form on C^" given by g = g o K7p^. It follows therefore that 

g(x,0 = Mx-e, (3.18) 

where M is a complex n x n matrix. We have 

H^=Mx-d,-M'^-d^, 

and using f ll.lOp . we see that the corresponding Hamilton map 

2 V -M7 

maps (x, 0) G Kt{A~^) to (l/2)(Mx, 0). Now the maps F and F are isospectral, and 
we conclude that, with the agreement of algebraic multiplicities, the following holds, 

Spec(M) = Spec(2F) n {Im A > 0}. (3.19) 

The final step in the normal form construction is the reduction of the matrix M in 
(I3.18P to its Jordan normal form. Such a reduction is implemented by considering 
a complex linear canonical transformation of the form 

Kc : 3 (x, ^ {C-'x, C'O e C'", (3.20) 

where C is a suitable invertible complex n x n matrix. On the operator level, 
associated to the transformation in fl3.20p . we have the operator 

Uc : u{x) H-> |det C\ u{Cx), (3.21) 

which maps the space iJ<i,o(C") unitarily onto the space iJ$^(C"), where $i(x) = 
^o{Cx) is a strictly plurisubharmonic quadratic form such that fi;c(A$y) = A$^. We 
notice that the strict convexity property 

<^>i{x)^\x\\ xgC", (3.22) 

remains valid. 

We obtain therefore the following result, which only is a slight reformulation of 
Proposition 2.1 of [15] and is closely related to the discussion in Section 3 of |25j . 
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Proposition 3.1 Let g : R " — ?> C be a quadratic form, such that Reg > 0. The 
quadratic differential operator 

Op"(g) : L2(R") ^ L2(R"), 

equipped with the domain 

ViOp'^iq)) = {ue L2(R"); {x' + {D,f) u G ^^(R")}, 
is unitarily equivalent to the quadratic operator, 

with the domain 

P(Op'^(g)) = {ue H^ACn-, (1 + e 4^(C")}. 

Here 

g(x,0 = Mx-e, 

where M is a complex n x n block-diagonal matrix, each block being a Jordan one. 
Furthermore, the eigenvalues of M are precisely those of2F in the upper half-plane, 
and the quadratic weight function $i(x) satisfies, 

^i{x)r^\x\\ xeC". 

The unitary equivalence between the operators Op^{q) and Op^{q) is realized by 
an operator of the form Ui o T o U2 : L^(R") -> Hq,^{CT'), where U2 is a unitary 
metaplectic operator on L^(R"'), while T is an FBI-Bargmann transform of the form 
(13.91) . (13.101) . and U\ is an operator of the form (I3.2ip . 

Let us recall now that Ai, . . . are the eigenvalues of the fundamental matrix F of 
q in the upper half-plane, repeated according to their multiplicity. It follows from 
Proposition 3.1 that the action of the operator Op'"(g) on Hq,^{C^) is given by 

n \ ^ n—1 

j=i j=i j=i 

Let 

$(x) = (1/2) xeC', 
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be the standard radial weight function. We observe that the weighted spaces 
(C") and iir$(C") contain a common dense subset, namely the space of holomor- 
phic polynomials on C". Indeed, it is well known that the normalized monomials 
form an orthonormal basis in if$(C"), while the density of holomorphic polyno- 
mials in if$^(C") has been explained in [15], and follows from the fact that the 
finite linear combinations of the generalized eigenfunctions of Op"'(g) are dense in 
L^(R"). When restricted to the space of holomorphic polynomials, the action of the 
operator Op"'(g) on Hq>^ (C") is clearly similar to the action of the closed densely de- 
fined operator Op"'(g) on the space H<s,{C^), the corresponding unbounded densely 
defined similarity transformation S : iy<i,^(C") — )■ i/<i,(C") being the identity map, 
with 'DlS) being the space of holomorphic polynomials on C". 

Assume now that the PT-symmetry condition, q{K{x), —k,*{^)) = q{x,^), {x,^) G 
R^", is maintained, and that the spectrum of the quadratic operator Op'"(g) on 
L^(R") is purely real. According to Remark 2.2, we then know that the eigenvalues 
Aj, 1 < J < n, are purely imaginary. If furthermore the nilpotent part in the Jordan 
decomposition of F vanishes, then (13.231) shows that 



is seen to be self-adjoint in the Bargmann space if$(C"), by means of a direct 
computation. This establishes the sufficiency part in Theorem II. 2 [ 

It only remains to verify the necessity in Theorem 11.21 To this end, assume that 
the PT-symmetric elliptic quadratic operator Op"'(g), Reg > 0, is similar to a 
self-adjoint operator, in the sense of the discussion in Section 1. It follows that for 
each A G Spec(Op'"(g)), the corresponding spectral subspace Ex C iS(R'^) consists 
entirely of eigenvectors. We conclude then from Proposition 3.1 that the nilpotent 
part in the Jordan decomposition of F vanishes. The proof of Theorem 11.21 is 
complete. 

Remark 3.3. If we drop the PT-symmetry assumption and merely assume that 
there are no Jordan blocks in the Jordan normal form of the fundamental matrix F, 
then it follows from (I3.24p that the operator Op"'(g), acting on if$(C'^), is normal, 
and therefore, the original quadratic operator Op'"(g) acting on L^(R'^) is similar 
to a normal operator. 




(3.24) 



which, when equipped with the domain. 



V{Op^{q)) = {u G H^iC^y, (1 + \xf)u G L|(C")}, 
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4 Example: quantum harmonic oscillator with qu- 
adratic complex interaction 

The purpose of this section is to illustrate Theorem II ■2l by applying it to the following 
two-dimensional quadratic Schrodinger operator, 

Op"'(g) = Dl^ + Dl^ + ujIxI + ujIxI + 2igxiX2. (4.1) 

Here Uj > 0, j = 1,2, ui 7^ ijJ2, and g E H. The operator Op'"(g) has been 
considered in the work [8] as a quantum model of a non-isotropic two-dimensional 
harmonic oscillator, perturbed by an additional quadratic interaction, with a purely 
imaginary coupling constant. 

The quadratic operator Op'"(g) is globally elliptic and PT-symmetric, with the 
corresponding involution given by k{xi,X2) = (— Xi,X2). In [8], the authors show 
that a certain method of separation of variables fails to work for Op"'(g) precisely 
when 

2g = ±iul-ul). (4.2) 

Furthermore, in the case when (14.21) holds, it is shown in [8] that the eigenfunctions 
of Op'"(g) do not form a complete set in L^(R^). 

By applying Theorem 1.2, here we shall show the following result. 
Proposition 4.1 The spectrum ofOp^{q) is real precisely when 

- \ujI - ujI\ < 2g < \ujI - ujI\ , (4.3) 
and Op"'(g) is similar to a self-adjoint operator if and only if 

- \ujI - ujI\ < 2g < \ujI - ujI\ . (4.4) 

Proof: The Weyl symbol q of the operator Op'"(g) is 

2 

q{x, = J2 fe' + ^1^?) + ^^9x1X2, (4.5) 
i=i 

and the corresponding fundamental matrix F is given by 

/ 1 \ 



2 V 



1 
-ul -ig 
\ -ig -ul 0/ 



(4.6) 
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A straightforward computation then shows that 



det (F - XI) = \^ + {ul + col)X^ + uj\ujI + 
with the four eigenvalues A being given by 

A2 = 



(4.7) 



(4.8) 



The expression under the square root sign is strictly less than ((w^ + bj\)l'£f'^ so 
that the eigenvalues are non-vanishing, as we aheady know since q is elliptic. It 
vanishes precisely when (14. 2 p holds, and from (14. 8 p we conclude that the eigenvalues 
are simple precisely when (14. 2 p does not hold. When (14. 2 p holds, then the spectrum 
of F consists of two double eigenvalues. 




Ai — A2 



1/2 



A.s — \a 



1/2 



(4.9) 



When (14. 3 p holds, we see that the eigenvalues A of F are on the imaginary axis, so 
thanks to (12. ip . we know that the spectrum of Op"'(g) is real. When (14. 3 p does not 
hold, the square root in (14. 8 p is non-real, so the eigenvalues A are off the imaginary 
axis, and Op^{q) has spectrum away from R, as was observed in Remark 2.2. 

From the discussion so far, we get the first statement in the proposition concerning 
the reality of the spectrum of Op"^(g), and also the fact that (14. 4p is a sufficient con- 
dition for Op"'(g) to be similar to a self-adjoint operator. It only remains therefore 
to show that Op'"(g) is not similar to a self-adjoint operator when 2g = ±(a;^ — 1^2)5 
and in view of Theorem 11.21 it suffices to show that Jordan blocks do occur in F 
then. In this case, we have two eigenvalues. 



A. 



±1 



1/2 



each of algebraic multiplicity two. If no Jordan block is present in the spectral 
subspace of A+, say, then the rank of F — A+J is equal to two. On the other hand, 
we have, using (14. 6p . 



F-\^I 



/-Ah 








-ig 

-ujI 



1 

1 

-A+ 

-A 



(4.10) 



+ 
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and the last three columns are linearly independent, since, 



det 












A 



+ 



1 


-A 



) 



The proof is complete. 



□ 



A Positivity 



The purpose of this appendix is to provide a brief but complete account of the 
relevant aspects of the theory of positive complex Lagrangian planes, required in the 
proof of Theorem II .21 See also Chapter 11 of [26], where a more general discussion 
is given. 

Let S C C^" be a real subspace of dimension 2n which is symplectic in the sense 
that (T|2 is a real and nondegenerate 2-form, where a = Yli d^j ^ dxj is the complex 
symplectic 2-form on C^"' = C" x C^. Then S is maximally totally real in the sense 
that S n z(Zl) = and the real dimension dimp{,(E) is maximal (2n). 

Let i = ly, '■ C^" — )■ C^" be the unique anti-linear map which is equal to the identity 
on E. Clearly, 



Example A.l. E = R^"', l{p) = p is the usual complex conjugation. 

Example A. 2. Let $ = be a real quadratic form on C" such that the Levi 
matrix d^dx^ is nondegenerate. Put 



Then S is symplectic. In fact, using x to parametrize E we get by a straightforward 
computation 



Using only that $ is real, we see that this restriction is real. Hence S is an I- 
Lagrangian manifold, i.e. a Lagrangian manifold for the real symplectic form Imcr. 



* — 

L a = a. 




(A.l) 




When the Levi form is nondegenerate, we see from the above expression that a\j, is 
nondegenerate. 
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Let be the unique holomorphic quadratic form on C^^y such that 

Differentiating this relation, we get 

keeping in mind that all second derivatives of $ are constant. 
The involution t = is given by 



i dy J \ i dx 



or more explicitly by 



Let A C C^*^ be a C-Lagrangian space, i.e. a complex n-dimensional subspace 
such that 0"!^ = 0. If S C C^" is a real symplectic subspace of maximal dimension 
with associated involution l = l\, we consider the Hermitian bilinear form 

b(p, fi) = -cr(p,L(fi)) on A X A. (A.3) 



The form ( ]A.3I) was introduced by L. Hormander [16], in the case when S = R 



2n 



Proposition A.3. The form ( \A.3\i is nondegenerate if and only if A (IT, = {0}. 

Proof: If 7^ p G A belongs to the radical of b, then l{p) is symplectically orthogonal 
to A, and consequently, G A. The vectors + and ^{p — l{p)) belong to 
S n A and at least one of them is 7^ 0, so S fl A 7^ {0}. Conversely S fl A is contained 
in the radical of b, so we have shown that the radical of b is non-zero precisely when 
SnA7^{0}. □ 

Example A. 4. If S = R^" and A is transversal to the fiber F = {(0,^; ^ G C"}, 
then A = A^: ^ = (f)'{x) = (f)"x, where is a complex holomorphic quadratic form 
on C". For {x,(f)'{x)) = (x, G A,^ we get 

^^(P.P) = ilra(p")x -x, 
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where 

Im0" = i(0"-(0"r), 
2i 

so the signature of the form (1A.3P is equal to that of Im0". 

Definition A. 5. Let A C C^" be a C-Lagrangian space and let S C C^" be a 
maximal real symplectic subspace with the corresponding involution l. We say that 
A is E-positive (or E- negative), if the Hermitian form (1A.3|) is positive (respectively 
negative) definite. 

Example A. 6. If E = R^" then A is E-positive if and only if A = A^, where 
Im0" > 0. 

Proof: If A = A(^ we already know by Example A. 4, that E-positivity of A is 
equivalent to that of Im0". This proves the if-part of the proposition. Conversely, 
if A is E-positive, then we see that A is transversal to the fiber F, so A = A^ and 
Example A. 4 applies again. □ 

The following result is a corollary of Proposition A. 3 and Example A. 4. 

Corollary A. 7 Let E he a fixed maximal real symplectic subspace ofC"^^. Then the 
set of all Tj-positive C-Lagrangian spaces is a connected component in the set of all 
C-Lagrangian spaces that are transversal to E. 

In fact, after applying a suitable complex canonical transformation, we may assume 
that E = R^"^ and we see from Example A. 6 that the set of all E-positive C- 
Lagrangian spaces is connected. Proposition A. 3 then shows that it is a connected 
component in the set of all C-Lagrangian spaces that are transversal to E. 

Now, let E = A$ be as in (lA.ip . where 

> 0. (A.4) 

Proposition A. 8. The fiber F = {(O,?]); r] E C"} is Tj-negative. 
Proof: Using (1A.2I) we see that ^(O,//) = {x,^) is given by 

i = X = (A.5) 



It follows that 



^a(0, r/; x, = -^r, ■ ($^,J-'r/ < ~\r^\\ 



□ 
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As we saw in (lA.Sp . the E-positive space l{F) is given by 



(A.6) 



Here, we notice that = $pih(a;) + ^hcimix), where 




($(x) - 



is the pluriharmonic (plh) part and 



-($(x) + ^{ix)) 



is the (positive definite) Hermitian part. We conclude that the positive C-Lagrangian 
space l{F) is of the form l{F) = A^^jj^, where $(a;) — ^p\h{x) x 

Proposition A. 9. Let S = A$ be as in flA.ip . f lA.4p . A C-Lagrangian space is 
Ti-positive if and only if A = Aj, where $ is pluriharmonic and $ — $ x 

Proof: If we decompose the pluriharmonic form $ as $ = $pih + $, we see that 
$ — $ > precisely when < ^hcrmix), x 7^ 0. Consequently, the set 

{Aj; $ is plh and $ — $ > 0} is a connected component of the set of all C- 
Lagrangian spaces that are transversal to Aij,. It contains A<j,pjj^ which is A<j,-positive, 
so by Corollary A. 7 it is equal to the set of all C-Lagrangian spaces, which are A$- 
positive. □ 

B Quadratic Fourier integral operators in the com- 
plex domain 

The Fourier integral operators encountered in this appendix arise when quantizing 
complex linear canonical transformations. The following discussion can therefore be 
viewed as a linear version of the theory presented in Chapters 3 and 4 in [26]. See 
also Chapter 3 of [20] . 

We start with the formal theory. A Fourier integral operator is an operator of 
the form 




(B.7) 
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where a G C and is a holomorphic quadratic form on C"^^"^^. We assume that 
is a nondegenerate phase function in the sense of [IT], so that 

d-^-— , d^— are hnearly independent, (Bi 



or equivalently that 



0',',(,,j^,,)isofrankiV. (B.9) 

We may assume without loss of generahty that 

ct>l, = Q. (B.IO) 

Indeed, if 7^ 0, we may assume after a hnear change of coordinates in 6', that 
e = {9', 9") e C^-*^ X C^, (p'^, Q = and that det 0'^'. g,, ^ 0. Then 9" ^ 
has a unique critical point 9"{x,y) which is nondegenerate and integrating out the 
^"-variables in (lB.7p . we get formally 

Au{x) = [[ e''^^^'y^^\u{y)dyd9\ (B.ll) 



where 6 is a constant non- vanishing multiple of a and ipi^x, y, 9') = (p{x, y, 9', 9'^) is a 
nondegenerate phase function with ipg, = 0. 

Assuming (IB.lOp until further notice, we see that (1B.9P becomes 

rank(0;;^(,_^)) = iV, (B.12) 

and in particular, < 2n. 
Let 

= {{x, y, e) e 0^"+^; 0;(x, y, 9) = 0} (B.13) 

This is a subspace of dimension 2n and we consider the corresponding canonical 
relation 

k: (y,r/)h^(x,0, (B.14) 

defined by its graph 

graph (k) = {(x, ^; y, r]) = (x, 0^.(a;, y, 9); y, -(f)'y{x, y, 9)); (x, y, 9) e C^}. (B.15) 
It is easy to check that dim (graph (k)) = 2n and that 
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where a^^^ — ay^^ = d^j A dxj — drjj A dyj. 
Assume that 

K is a canonical transformation, (B.16) 

which means that the maps graph (k) 3 {x,^]y,ri) G C^" and graph (k) 3 

{x,C,] y, rj) {y, r]) G C^" are bijective. Actually, the bijectivity of one of the maps 
implies that of the other. In fact, the bijectivity of the first map is equivalent to 
that of 3 (x, y, 9) I— J- (x, y,9)) G C^" which is equivalent to that of 

0^"+^ 3 (x, y, 9) ^ (x, 01., = (x, ,.x + <Pl^y + ct^%9, ct^^^x + 0;',^|/) G 0^"+^. 

This is equivalent to the bijectivity of 

and implies that 0g ^ and 0" g are of maximal rank N and in particular that N < n. 
We can write 

N 

0(x, y, 9) = g{x, y) + J2 y)^v (B-IT) 

i=i 

where is quadratic and fj are linear. Further, d^fi, ....d^fN are linearly indepen- 
dent and similarly for the y-differentials. The subspace is given by the equations 
fj{x,y) = 0, j = 1, .., in C^'^"'"^ and the same equations in the x,y space define 
the projection of C<^, or equivalently, the projection tt^j^j^ (graph (k)), so A^ is uniquely 
determined by k. Let 

N 

0(x, y, 9) = gix, y) + U^^ y)^J (B-18) 

be a second nondegenerate phase function that generates the same canonical trans- 
formation K. Then /i = ... = /at = and^/i = ... = /jy = define the same subspace 
of C^" and after a linear change of the 9 coordinates, we may assume that 

N 

(f) = (f){x,y,9) = g{x,y) + ^fj{x,y)9j. 

i=i 

The fact that the two phases give rise to the same canonical transformation now 
means that there are linear forms kj{x, y), I < j < N such that 

N 

9{x, y) = g{x, y) + ^ fj{x, y)kj{x, y). (B.19) 
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Then 

(j){x,y,e) = <j){x,y,e + k{x,y)), (B.20) 

which will imply that (p gives rise to the same Fourier integral operators, once we 
have explained how to choose the contour of integration in (IB .71) . 

Let $0 a-iid $i be real quadratic forms on (different copies of) C". 

Proposition B.l Let (p be a nondegenerate phase function as above, with the as- 
sociated canonical transformation n. The following statements are equivalent: 

1) k(A$J = A$^. 

2) The quadratic form 

{y,e)^%{y)-\m<P{^,y,e) (B.21) 

is nondegenerate so that 

{y,e)^%{y)-\m<P{x,y,e) (B.22) 
has a unique critical point {y{x),6{x)) for every x G C". Moreover, 
$i(x) = vcj,,0($o(y) -Im0(x,?/,^)). 

Proof: We first assume 2). Then at the critical point {y{x),9{x)), we have 
^<f){x,y,e) = 0, 



09 

2a$i 2d 



{-lm(j)){x,y,9) = — (x,y, 



i dx i dx dx 

which means that 1) holds. 

Now assume 1). Then for every x G C" there exists a unique y G C" such that 
X G Trx{K,{y, l^^(x))). Equivalently for every x G C", ( IB. 221) has a unique critical 
point {y{x),6{x)). Since we are dealing with second order polynomials this critical 
point is nondegenerate. If $i(x) is the critical value, we also see that 

2 _ 2 d^i 

i dx i dx 
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and since $i and $i are quadratic forms, we conclude that $i = $i. □ 

If k(A$o) = A<i,j and $0 is pluriharmonic (so that A$q is a C-Lagrangian space), 
then $1 is also pluriharmonic and the nondegenerate quadratic form flB.2ip has the 
signature {n + N ,n + N). In general, when $0(1/) and hence $0(2/) — Ini0(O,?/,6') 
is plurisubharmonic, then flB.2ip has at most n + N negative eigenvalues, and in 
order to find a suitable integration contour, when realizing A, we wish this maximal 
number of negative eigenvalues to be attained. The following result gives a positive 
answer in the case that we are interested in. 

Proposition B.2 Assume that K(A<i,y) = A<j,^, where d^dz^Q and d^dz^i are posi- 
tive definite. Then the signature of the quadratic form flB.211) is equal to + A^, n + 
N). 

Proof: Let p be a positive definite quadratic form on A$q and let p also denote the 
holomorphic extension to C^" (which is a well defined holomorphic quadratic form 
since A$g is maximally totally real). Then the spectrum of the fundamental matrix 
Fp is given by the eigenvalues ±iyUi, ±ifin, where fij > 0. Let Ao,+oo denote the 
spectral subspace of C^" corresponding to ifii, Then we know from p5] that 
Ao,+oo is positive with respect to A^^, so that, according to Proposition A. 9, Ao,+oo 
is of the form 1^ = ^^^q^{x), where $0,4-00 is a pluriharmonic quadratic form such 
that ^o{x) - ^o,+oo{x) X Ixp. 

The spectral subspace Ao,+oo can be viewed as the unstable manifold for the H^ip- 
flow and if we put Aq,* = exp (— 'ji:ifp)(A$g), for t G R, then Aq,* is I-Lagrangian 
and R-symplectic, and Aq,* converges to Ao,+oo exponentially fast. Furthermore, at 
least for small t > 0, we have Ao,f = A<j,g^, where $o,t(2^) is the real quadratic form 
obtained by solving the eikonal equation 

^$o,t(x) + Rep (^x, -^^<i>o,t{x)^ = 0. (B.23) 

See for instance, [27], as well as Remark 11.7 in [13], for more details. Notice here 
that since p is constant along the H^ip trajectories, we have Rep = p in (IB. 231) and 
P (•^' f ^^0'*("^)) ^ I'^l^' iion-uniformly with respect to t. 

A priori it is not clear that we can solve f lB.23p for all t > 0. However, as long 
as the solution exists, we see that t H- $o,t is a decreasing family of strictly 
plurisubharmonic quadratic forms, so that $o,t(a^) < ^o{x). For such a t, we write 
Ip := $o,t = V^herm + V'pih, wherc ^pherm{x) = l{ip{x) + ip{ix)) IS the Hermitian 
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part and ippi^ix) = ^{ip{x) — ip{ix)) is the pluriharmonic part. Similarly, we write 

:= $o(x) = 0herm + 0plh- FrOIIl ijj < (j) WS get < Iphcvm < li<Pix) + = 0herm- 

Also, 

ll^plhix) < (0herm(x) - 1phcrm{x)) + (f)pll,{x) < 

-^pih(a;) = ^pih(^a;) < 0(zx), 

so that 

< V^pihlx) < (p{x). 

This means that we have t-independent upper and lower bounds on $o,t) which 
prevent an explosion in the eikonal equation flB.23p . Consequently, Ao,t = A$g^ for 
all < t < oo, and the exponential convergence of Ao,t to Ao,+oo when t — )■ +00 
shows that 

$o,t '^'o.+oo exponentially fast, when t — )■ +00. 

Let us define q on C^" hj q o k, = p. Replacing ($05^) in the above discussion by 
($1, q), we get a decreasing family of strictly plurisubharmonic quadratic forms 
for < t < +00, converging to the pluriharmonic form $i,+oo when t — > +00, with 
the property that 

fi;(A$Q_J = A<i,^^, < t < +00. 
Proposition B.l now shows that 

is a nondegenerate quadratic form for < t < +00, necessarily with t-independent 
signature. When t = +00 the signature is equal to {n + N,n + N) so this is also the 
case when t = 0, which gives the proposition. □ 

Now let us recall the representation (lB.17p . where (p'^ y and (p^ .^ are of maximal rank 
A^. After separate linear changes of the x and y coordinates, we may assume, in 
order to simplify the notation only, that 

N 

<P{x, y, 9) = g{x, y) + ^ {x, - yj)e, = g{x, y) + [x" - y") ■ 6, (B.24) 

j=n-N+l 

writing x = {x', x") G C"~^ x and similarly for y. We see directly that 
-Im 0(0, 0, y", 9) + $0(0, y") = -Im (7(0, 0, y") + %{0, y") + Im {y" ■ 9) 
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is a nondegenerate quadratic form of signature {N,N). More generally, the second 
order polynomial 

iy", 6) ^ -Im y, 6) + $0(1/) = -lmg{x, y) + <l>o(y) - Im ((x" - y") ■ 6) (B.25) 

has a unique critical point of signature (2A^, 2A^). The critical point {y'^{x, y'), Oc{x, y')) 
is given by y'^. = x", 6c = j^friy' , x") + §^{x, y' , x") and the critical value is 

-lm(t){x,y\x'\e,) + ^^{y\x") = -\mg{x,y' ,x") + ^o{y\x"). (B.26) 

On the other hand, by Proposition B.2, we know that 

(y,^)^-lm0(x,y,^) + $o(y) (B.27) 

has a unique critical point which is of signature {n + N,n + N). It follows that 

y' ^ -\mg{x,y\x") + ^Q{y\x") (B.28) 

has a unique critical point which is of of signature {n — N,n — N), and the corre- 
sponding critical value is equal to that of (lB.27p . namely $i(x). 

Let 7 C Cy~^ be a real-linear subspace of dimension n — N such that 

-Im(7(0,2/',0) + <l>o(y',0)x-|yf, G 7. (B.29) 

Then 

- \mg{x,y\x") + $o(y',a:") - $1(0:) x -\y' - y[{x)\\ y' E + 7. (B.30) 
Let us consider the contour r(x) C C^^^ of real dimension n + N, given by 

y' G y'cix) + 7, y" G C^, ^ = y') + zC{x"-y"), 
parametrized by {y',y") G 7 x C^. Along T{x) we have 

- Im0(x, y, 6) + $0(2/) - $i(x) x - - y'^{x)\^ + |/ - a:'f ) , (B.31) 
provided that we choose C > large enough. 

Let u G if<i,„(C"), and let us realize Au = Aj-u in ( IB. 71) by integrating over T{x). 
We get 

Jy"£C^ J y' &y'^{x)+-i (B.32) 

= e-*Mre*«(Me-*«), 
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where — Im ipri^, u) is equal to the right hand side in (lB.3ip . From (lB.3ip we see that 
the integral in (1B.32P converges for each x G C", and by a standard application of 
Stokes' formula, we also know that if F is a second contour with the same properties 
as F, then Ayu{x) = Aj;u{x). 

We now wish to estimate the C^L"^, L^)-norm of the effective operator e~*^ Ape*" in 
(lB.32p . We know that the map 

X ^ yc{x) = {y'cix), x") (B.33) 

is bijective. Replacing x by the new x coordinates {y'^{x) , x") , we may reduce our- 
selves to the case when y'^lx) = x'. After a simultaneous rotation in the x' and y' 
variables we may further assume that 7 = R*^"^. Then from flB.3ip . flB.32p . we get 

(B.34) 

By Schur's lemma, 

||e-*Mrn|U.(^,)<0(l)||e-*"n|U2(^^,), 

uniformly for t' G R*^"^, where Ef = {x G C"; Imx' = t'}. Integrating this with 
respect to t', we get 

||e-*MrM|U2(c.) < 0(l)||e-*«M|U2(cn). 
We have proved the following result. 

Proposition B.3 Realizing the operator A by means of the contour T above, we get 
a bounded operator A = Ar : i?$o — )■ . We get the same operator if we replace 
F by any other contour with the same properties. 



The next result can be proved by the general methods of |26j . 

Proposition B.4 Let (-B,k) have same general properties as {A, n) and assume 
that k(A$J = A$2, where $2 is a strictly plurisubharmonic quadratic form. Let F 
be a contour allowing to realize B = : H<s>^ — t- Hi^.^. Then there exists a Fourier 
integral operator C as above with associated canonical transformation H o n and a 
contour F such that Cf = Bf o A-p. 
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Proof: We shall only recall the general ideas of the proof. First of all, the elimination 
of superfluous 6 coordinates above was only in order to simplify the presentation, 
and it is very easy to realize Fourier integral operators with more than the minimal 
number of such fiber variables. Formally Bo A then becomes such a Fourier integral 
operator with k o k as the associated canonical transformation. When composing 
Bf o Ar we automatically get a good contour, and as we have pointed out, this 
contour can be deformed into any other good contour for B o A. Superfluous 6 
variables can then be eliminated by the exact version of stationary phase. □ 

A special case is when k = id. We then get a multiple of the identity operator. If k 
in the last proposition is equal to then we conclude that B o A is a. multiple of 
the identity, and assuming that the amplitude a in (]B.7p is 7^ 0, we can choose the 
amplitude 6 in a corresponding representation of B in such a way that B o A = I . 
In particular, we see that Ap : if^^ — > in Proposition B.3 has a bounded 
(two-sided) inverse. 
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